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INTRODUCTION 

In  the  last  two  decades  the  methods  outlined  by  St.  Denis  and 

1* 

Pierson  for  the  estimation  of  the  magnitude  of  oscillatory  ship  motions 
in  irregular  seas  have  become  firmly  established  in  engineering  practice. 
These  methods  apply  strictly  only  to  ship  responses  which  can  be  assumed 
to  be  a linear  function  of  wave  height,  and  they  involve  as  well  the  as- 
sumption that  the  wave  process  is  Gaussian. 

In  this  context  the  first  recognized  ship  motions  problem  (ship 

rolling)  remains  a problem  to  some  extent.  Roll  damping,  at  least  for 

low  or  zero  ship  speeds,  has  been  considered  to  be  a mixture  of  linear  and 

2* 

quadratic  damping  for  about  a century.  No  modern  hydrodynamic  analysis 
has  challenged  the  model  — for  that  matter  there  appears  to  be  no  theoret- 
ically based  prediction  method  for  roll  wnich  is  completely  free  of  empiri- 
cism with  respect  to  roll  damping.  Modern  theory  appears  to  consistently 
underestimate  the  linear  part  of  roll  damping  for  low  speeds.  However,  if 
the  necessity  for  empiricism  is  accepted  (as  it  generally  is,  and  will  be 
herein)  there  is  still  another  problem  with  roll  and  this  is  the  simple 
fact  that  a nonlinearity  often  appears  to  exist  for  low  to  moderate  rolling 
ampl i tudes. 

In  the  majority  of  applications  to  design  what  is  most  wanted  is 
a measure  of  the  statistics  of  rolling  maxima  to  be  expected  in  various 
hypothetical  or  realistic  irregular  sea  conditions.  When  frequency  domain 
predictions  of  the  statistics  of  roll  in  irregular  seas  are  made,  especially 
in  a multi-degree  of  freedom  problem,  the  usual  approach  is  to  linearize 
the  damping  coefficient  (compensating  at  the  same  time  for  the  uncertainties 


* 

1.  St.  Denis,  M.  and  Pierson,  W.J.,  Jr.,  On  the  Motions  of  Ships  in 
Confused  Seas,"  SNAME  Vol.  6l,  1953. 

2.  "The  Papers  of  William  Froude,  "The  Institution  of  Naval  Architects, 
London , 1955* 
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involved  in  the  estimation  of  the  linear  part)  and  thereafter  to  conduct 
the  analysis  and  make  the  predictions  as  though  the  system  were  completely 
linear.  Alternately,  rolling  non-linearities  can  be  (and  have  been) 
incorporated  into  time  domain  solutions  for  roll  response  to  irregular 
waves.  In  this  case  the  validity  of  the  statistical  prediction  of  maxima 
depends  to  a great  extent  upon  the  amount  of  sample  of  computed  response 
which  is  generated  --  and  consequently  upon  the  time  and  money  available 
for  the  solution  to  a particular  problem.  Accordingly,  in  most  cases 
involving  the  linear  frequency  domain  approach  the  computation  process  is 
economical  but  the  results  may  have  to  be  interpreted  with  care  because 
a basic  assumption  of  the  prediction  framework'  is  violated.  In  contrast, 
in  the  time  domain  approach  there  are  fewer  problems  about  the  adequacy 
of  representation  but  the  costs  tend  to  be  considerably  greater  for  equiva- 
lent results. 

Ideally,  what  would  be  very  useful  in  the  prediction  of  ship  rolling 
statistics  in  random  seas  is  a prediction  framework  analogous  to  the  linear 
framework  of  St.  Denis  and  Pierson'  in  which:  a)  at  least  weak  non-linearities 

could  be  accommodated  for  the  general  multi-degree  of  freedom  situation, 
b)  multi-directional  seas  could  be  considered  as  input,  c)  the  hydromechanic 
data  required  could  be  produced  with  conventional  techniques,  d)  the  pre- 
diction or  a major  part  could  be  carried  out  in  the  frequency  domain  for 
economy  (as  well  as  to  take  advantage  of  the  accumulating  frequency  domain 
descriptions  of  real  sea  waves)  and  e)  the  statistics  of  maxima  could  be 
estimated  with  firmly  based  theory  in  which  the  possible  effects  of  non- 
linearities  are  accounted  for. 

It  is  clear  that  this  goal  is  not  in  hand,  and  equally  clear  that 
if  the  required  technology  exists  it  is  unlikely  to  be  reduced  to  practice 
in  the  immediate  future.  However,  a number  of  steps  in  this  general 
direction  have  been  made,  and  it  was  the  general  objective  of  the  present 
work  to  add  another. 

Previous  work  on  the  problem  of  a frequency  domain  prediction  of 
non-linear  rolling  in  irregular  seas  has  centered  upon  the  problem  of 
predicting  the  spectrum  and  the  variance  of  zero  speed  rolling,  this  being 
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the  case  in  which  the  damping  non-linearity  has  been  found  to  be  most 
obvious,  and  the  case  for  which  it  is  plausible  to  reduce  the  problem  to 

a single  degree  of  freedom.  Typically,  a single  degree  of  freedom  rolling 

2 

equation  not  far  different  from  that  of  W.  Froude  is  solved  in  some  sense 

3*  4* 

for  the  random  excitation  case.  The  work  of  Kaplan  and  Vass i lopoulos 
involve  equivalent  linearization  techniques  in  the  estimation  of  the 
variance  of  roll  and  the  roll  spectrum.  In  this  technique  nonlinear 
elements  are  replaced  by  linear  elements  chosen  so  as  to  minimize  the 
resulting  mean  square  errors  for  the  case  of  random  excitation.  Yamanouchi^ 
approached  the  problem  with  a perturbation  technique  for  the  solution  to 
the  non-linear  differential  equation.  In  this  solution  the  spectrum  of 
non-linear  rolling  turns  out  to  be  the  sum  of  the  linear  roil  spectrum 
and  various  convolutions  of  the  linear  velocity  spectrum. 

For  practical  purposes  the  result  of  both  the  above  approaches  is 
the  predicted  roll  variance.  The  implicit  assumption  is  that  the  statistics 
of  maxima  are  adequately  described  by  the  Rayleigh  distribution  with  para- 
meter equal  to  square  root  of  variance.  This  assumption  has  been  partially 
vindicated  by  a study  in  which  numerical  simulation  of  the  single  degree 
of  freedom  rolling  equation  was  carried  out.  The  results  indicated  that 
while  the  Rayleigh  distribution  is  probably  not  a completely  proper  as- 
sumption, reasonably  good  predictions  under  this  assumption  could  be 
expected  for  quantile  averages  up  to  average  of  1/10  highest  amplitudes, 
despite  inclusion  of  non-linearities  within  the  nominal  range  of  magnitude 
observed  in  unstabiiized  ships  and  models. 


3.  Kaplan,  P.,  "Lecture  Notes  on  Non-Linear  Theory  of  Ship  Roll  Motion 
in  a Random  Sea  Way, " ITTC  Transactions,  1966. 

* 

4.  Vass i lopoulos,  L.,  "Ship  Rolling  at  Zero  Speed  in  Random  Beam  Seas 
with  Non-Linear  Damping  and  Restoration,"  Journal  of  Ship  Research, 
Vol.  15,  No.  b,  December  I97I. 

* 

5.  Yamanouchi,  Y.,  "On  the  Effects  of  Non-Linearity  of  Response  on  Cal- 
culation of  the  Spectrum,"  ITT C Transactions,  1966. 

* 

6.  Dalzell,  J.F.,  "A  Note  on  the  Distribution  of  Maxima  of  Ship  Rolling," 
Journal  of  Ship  Research,  Vol.  17>  No.  4,  December  1973. 
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A third  basic  approach  to  the  prediction  of  the  statistics  of  non- 

7* 

linear  roiling  is  afforded  by  the  Fokker-Planck  equation  method,  Caughey  . 
Evaluation  of  roll  response  stat i st i cs  accord i ng  to  this  approach  promise 
to  be  extremely  difficult  when  the  spectrum  of  excitation  is  not  white 
(flat).  Haddara  uses  a modification  of  this  approach  resulting  in  esti- 
mates of  roll  variance  for  the  case  when  the  excitation  is  white. 

A fourth  basic  approach,  and  the  subject  of  the  present  work,  is 
the  functional  series  model.  The  general  attractions  of  the  approach  are 
several.  Among  these  are  that  as  a conceptual  framework  the  model  is 
suitable  for  any  reasonably  well  behaved  wave  input  (regular,  transient, 
or  random),  and  since  it  contains  the  completely  linear  system  as  a 
special  case  it  appears  to  have  the  potential  of  being  a logical  extension 
to  present  practice.  In  addition,  theoretical  prediction  methods  for 
spectra  may  be  derived,  and  it  appears  that  it  may  be  possible  to  approxi- 
mate the  statistics  of  maxima.  Finally,  it  is  possible  in  principle  to 
relate  the  functions  required  by  the  model  to  the  results  of  hydromechani- 
cal analyses  and  experiment. 

Wiener  introduced  the  functional  series  model  (also  called  Volterra 

Series)  into  non-linear  circuit  analysis  during  World  War  II,  and  sub- 

9* 

sequently  published  this  work  and  extensions  . Over  the  intervening 
years  the  ideas  and  applications  have  slowly  been  amplified  and  simplified 
for  consumption  in  the  communication  and  electronics  fields.  Of  the  many 


7.  Caughey,  T.K.,  "Derivation  and  Application  of  the  Fokker-Planck  Equa- 
tion to  Discrete  Nonlinear  Dynamic  Systems  Subjected  to  White  Noise 
Random  Excitation,"  Journal  of  the  Acoustical  Society  of  America, 

Vol.  35,  No.  11,  November  1963. 

* 

o.  Haddara,  M.R.,  'A  Modified  Approach  for  the  Application  of  Fokker- 
Planck  Equation  to  the  Nonlinear  Ship  Motions  in  Random  Waves," 
International  Shipbuilding  Progress,  Vol.  21,  No.  2142,  October  197^. 

* 

9.  Wiener,  N.,  Nonlinear  Problems  in  Random  Theory,"  The  Technology 
Press  of  MIT  and  John  Wiley  and  Sons,  Inc.,  1958. 
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. 10*  11* 
papers  in  that  literature,  those  of  Barrett  and  Bedrosian  and  Rice 

may  be  recommended. 

Application  to  seakeeping  problems  was  suggested  indirectly  by 

12*  . 13*  . )ii* 

Tick  in  I96I.  Hasselman  and  Vassi lopoulos  indicated  direct 

applications  to  some  classes  of  seakeeping  problems  roughly  10  years  ago. 

. . 1 5*  16*. 17* 

More  recently  Neal  and  Dalzell  ’ have  presented  results  bearing 
on  the  application  of  the  functional  series  model  to  what  is  probably  the 
simplest  non-linear  seakeeping  problem,  the  ship  resistance  added  by  waves. 

In  the  application  to  ship  rolling  there  is  a serious  obstacle. 

14 

This  was  pointed  out  by  Vassi lopoulos  . It  is  that  if  the  functional 
expansion  and  a differential  equation  are  to  be  related,  it  appears  that 
the  equation  must  be  analytic  for  small  values  of  the  variables.  This  is 
not  the  case  for  the  "quadratic"  term  ordinarily  used  to  represent  the 
damping  non-linearity. 

•X- 

10.  Barrett,  J.F.,  "The  Use  of  Functionals  in  the  Analysis  of  Non-Linear 
Physical  Systems,"  Journal  of  Electronics  and  Control,  Vol-  I 5,  No.  6, 
December  1963. 

11.  Bedrosian,  E.  and  Rice,  S.O.,  "The  Output  Properties  of  Volterra 
Systems  (Nonlinear  Systems  with  Memory)  Driven  by  Harmonic  and  Gaussian 
Inputs,"  Proceedings  of  the  IEEE,  Vol.  59»  No.  12,  December  1971* 

* 

12.  Tick,  L.J.,  "The  Estimation  of  the  'Transfer  Functions'  of  Quadratic 
Systems,"  Technometrics,  Vol.  3,  No.  4,  I96I. 

* 

13-  Hasselman,  K. , "On  Non-Linear  Ship  Motions  in  Irregular  Waves," 

Journal  of  Ship  Research,  Vol.  10,  No.  1,  1966. 

* 

14.  Vass i lopoulos,  L.A.,  "The  Application  of  Statistical  Theory  of  Non- 
linear Systems  to  Ship  Motion  Performance  in  Random  Seas,"  Inter- 
national Shipbuilding  Progress,  Vol.  14,  No.  150,  I967. 

* 

15*  Neal,  E.,  "Second-Order  Hydrodynamic  Forces  Due  to  Stochastic  Excita- 
tion," Tenth  ONR  Symposium  on  Navai  Hydrodynamics,"  Massachusetts 
Institute  of  Technology,  197^- 

* r 

16.  Dalzell,  J.F.,  "Cross-Bispectral  Analysis:  Application  to  Ship 

Resistance  in  Waves,"  Journal  of  Ship  Research,  Vol.  18,  No.  1, 

March  1974,  pp  62-72. 

* 

17.  Dalzell,  J.F.,  "Application  of  the  Functional  Polynomial  Model  to 
the  Ship  Added  Resistance  Problem,"  Eleventh  Symposium  on  Naval 
Hydrodynamics,  University  College,  London,  1976. 
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However,  because  of  reasonable  success  achieved  with  the  appli- 
cation of  the  functional  expansion  method  to  the  added  ship  resistance 

problem,  the  possibilities  of  application  to  ship  rolling  were  re- 

. 18* 
investigated.  This  report  and  a companion  report  are  the  results 

of  the  re-investigation.  Relative  to  the  ideal  goals  previously  cited, 

the  objectives  of  the  study  were  much  simplified.  Briefly,  the  study 

objectives  were  limited  so  as  to  correspond  with  the  objectives  and 

restrictions  of  previous  work^,^,^,^j  that  is,  the  objective  of  the 

present  work  was  to  attempt  to  apply  the  functional  series  model  to  the 

estimation  of  the  spectrum  of  roll  response  to  random  excitation  under 

the  assumption  that  the  usual  single  degree  of  freedom  roll  equation  holds. 
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18.  Dalzell,  J.F.,  "A  Note  on  the  Form  of  Ship  Roll  Damping," 
S I T-DL-76- 1887,  Davidson  Laboratory,  Stevens  Institute  of 
Technology,  May  1976. 
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THE  ASSUMED  SINGLE  DEGREE  OF  FREEDOM  ROLLING  EQUATION 

For  purposes  of  initial  discussion,  an  equation  taken  to  represent 
zero  speed  ship  rolling  in  beam  waves  may  be  written: 

1*9  + N(  cp)  + B(cp)  = F(t)  (1) 

where:  cp  = roll  angle 

I = rol 1 inertia 
N(cp)  = a roll  damping  function 
B(cp)  = a roll  restoring  function 
F(t)  = an  excitation  function  of  time  (t) 

Within  the  range  of  analytical  representations  of  inertia,  restoring 
moment  and  exciting  moment  thus  far  seen  in  conjunction  with  single  degree 
of  freedom  rolling,  no  fundamental  problems  of  application  of  the  function- 
al series  were  envisioned.  As  noted  in  the  introduction,  it  is  the  usual 
representation  of  the  damping  function  which  is  the  problem,  and  this  may 
be  written  as: 

N(  <p)  = N cp  + N Icplcp  (2) 

21  2 2 

Before  anything  much  is  worthwhile  in  applying  the  functional  series  to 
ship  rolling,  the  quadratic  term  must  be  dealt  with.  Because  this  term 
is  odd  in  cp,  the  simplest  approach  appeared  to  be  to  replace  the  term  with 
an  odd  series  in  cp.  This  possibility  was  investigated  in  some  detail  and 
the  results  appear  in  Ref.  18. 

The  approach,  results  and  conclusions  of  Ref.  18  may  be  summarized 
as  follows.  Physically,  most  of  the  explanations  of  the  origins  of  the 
form  (Eq.  2)  for  damping  moment  more  or  less  follow  W.  Froude's  views 
that  the  linear  term  is  caused  by  energy  dissipation  in  waves  made  by  the 
ship,  and  that  the  quadratic  term  represents  the  real  fluid  effects. 
Theoretical  estimation  methods  are  lacking  for  the  most  part.  Accordingly, 
nearly  the  entire  justification  for  the  linear-plus-quadratic  representa- 
tion for  roll  damping  is  that  it  appears  to  work  — in  the  sense  that  curves 


R-l8$)4 


of  declining  angles  obtained  in  ship  or  model  experiments  may  be  reasonably 
well  fitted  under  this  assumption.  Realistic  estimates  of  N21  and  N22  in 
Eq.  (2)  are  almost  totally  empirical,  and  those  which  exist  for  particular 
ships  come  overwhelmingly  from  sallying  experiments. 

The  general  approach  in  Ref.  18  was  therefore  also  totally  empiri- 
cal and  it  consisted  of  reanalyzing  a number  of  sets  of  ship  and  model 
sallying  results  under  the  assumption  that  the  damping  function  of  Eq.  (l) 
could  be  represented  in  the  form: 

. . .3 

N(cp)  = N cp  + N cp  (3) 

31  33 

In  this  effort  the  same  additional  assumptions  were  made  as  are  usual  in 
analyses  of  sallying  data  according  to  the  quadratic  model,  Eq.  (2).  These 
are  that  the  excitation  in  Eq.  (l)  is  zero,  the  inertia  is  a constant,  and 
that  the  restoring  function  is  linear  and  equal  to  AGMcp,  where  A is  ship 
displacement  and  "GM  is  transverse  metacentric  height.  The  results  of  this 
effort  were  that  within  the  range  and  probable  precision  of  experimental 
data  there  is  not  much  to  choose  between  the  two  damping  representations 
(Eqs.  2 and  3).  Outside  the  range  of  experimental  data  both  are  extra- 
polations, either  may  or  may  not  be  reasonable.  Given  the  possibility 
of  selecting  the  coefficients  for  each  model  in  a "best"  way,  either  model 
appears  capable  of  realistically  representing  experimental  curves  of  de- 
clining angles  or  the  curves  of  roll  extinction  derived  therefrom. 

Two  approaches  to  the  problem  with  the  quadratic  term  in  Eq.  (2) 
appear  reasonable.  One  is  to  use  Eq.  (3)  as  the  damping  function  under 
the  assumption  that  empirical  values  for  N31  and  N33  are  available  or  can 
be  developed. 

The  other  approach  developed  in  Ref.  18  is  to  replace  Eq.  (2)  by 
the  following  approximation: 
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This  approximation  results  from  a least  square  fit  of  a two-term  odd 
series  in  cp  to  |cp|cpover  the  roll  velocity  range  (-R  < cp  < R).  Accordingly, 
the  approximation,  Eq.  (4),  may  be  thought  of  as  having  been  designed  for 
roll  velocities  up  to  (R).  The  work  of  Ref.  18  implies  two  possible  points 
of  view  in  applying  this  approximation.  The  first  involves  fixing  R equal 
to  the  highest  roll  velocity  involved  in  the  experiments  from  which  the 
values  of  N21  and  N22  were  derived.  Under  these  circumstances  the  devia- 
tions of  Eq.  (4)  from  the  quadratic  form,  Eq.  (2),  for  1 cpl  s R are  of  the 
same  order  as  experimental  scatter  and  the  maximum  deviation  can  be  ex- 
pected to  be  2 or  3$  of  N(R),  for  typical  mixes  of  linear  and  quadratic 
damp i ng . 

The  other  point  of  view  involves  considering  the  quadratic  repre- 
sentation as  absolute,  and  its  coefficients  exact.  In  this  case  the 
approximation,  Eq.  (4)  can  be  thought  of  as  an  equivalent  non-linearization 
for  some  selected  roll  velocity  range,  ± R.  In  the  sense  of  making  esti- 
mates of  roll  variance,  this  implies  an  iterative  scheme  whereby  the 
approximate  range  of  roll  velocity  would  be  computed  for  an  assumed  value 
of  R and  the  procedure  repeated  as  required  to  insure  that  the  computed 
and  assumed  velocity  ranges  were  approximately  the  same.  Under  these 
conditions  the  maximum  deviations  from  the  quadratic  form,  Eq.  (2)  can 
also  be  expected  to  be  2 or  3$  of  N(R)  or  less  for  typical  mixes  of  linear 
and  quadratic  damping. 

To  put  some  perspective  upon  the  percentages  just  quoted  it  may  be 
pointed  out  that  the  approximation  to  the  quadrat ic  part  of  Eq.  (2)  is 
always  within  4%  of  N22RS,  the  2 or  3$  figures  quoted  for  N(  cp)  as  a whole 
result  from  consideration  of  typical  experimental  data.  If  the  quadratic 
part  of  Eq.  (2)  is  1 i near i zed  in  a least  square  sense  over  the  range 
( -R  < cp  < R),  the  quadratic  term  in  Eq.  (2)  would  be  replaced  by 
(O.T5N32R43)>  and  thus  a 1 inearizat  ion  of  the  quadratic  term  involves 
deviations  from  the  original  of  up  to  25$  of  N22R2.  Thus  even  when  the 
position  is  taken  that  the  quadratic  form,  Eq.  (2),  is  cast  in  concrete, 
the  "equivalent  non- 1 inear i zat ion"  results  in  a six-fold  improvement  over 
linearization  in  the  representation  of  the  non-linear  damping  term.  The 
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differences  are  shown  graphically  in  Figure  l for  positive  9.  Within  the 
range  of  roll  velocities  bounded  by  ±R  it  seems  clear  that  the  approxima- 
tion of  Eq.  (4)  represents  a significant  improvement  over  simple  lineari- 
zation. 

No  matter  which  of  the  aforementioned  approx imat ions  to  the  damping 
function  are  employed, the  form  is  the  same,  so  that  for  the  purposes  of 
proceeding  further  the  damping  function  will  be  assumed  to  be: 

.3 

N(cp)  = I o acp  + [a  9 /a  (5) 

1 3 


where: 


a2  = AGM/ I , the  squared  undamped  roll  frequency 


°i  = N3i  / 1 o 
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The  N.  . are  empirical  as  defined  in  Equations  (2)  and  (3),  and  the  alter- 
native forms  for  the  Cr^  correspond  to  the  various  approximations  just 
discussed. 

The  coefficient  of  9 in  Eq.  (l),  the  inertia,  is  conventionally 

assumed  to  be  a constant.  In  some  of  the  previous  work^’^  the  restoring 

function,  B(9)  has  been  assumed  to  be  linear  as  was  done  by  Froude.  In 

4 6 8 14  19* 

other  work  ’ ’ ’ * , the  restoring  function  has  been  taken  as  a two- 

term  odd  series  in  9,  that  is: 


J 
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19.  Haddara,  M.R.,  "On  Nonlinear  Rolling  of  Ships  in  Random  Seas," 
International  Shipbuilding  Progress,  Vol.  20,  No.  230,  October 
1973. 
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B(  cp)  = AGM(  cp  + Cep  ) 

Since  this  expression  is  of  the  same  order  as  Eq.  (5),  additional  terms 

6 

do  not  seem  justified,  though  an  additional  term  proportional  to  cp  would 
allow  a quite  good  representation  of  the  static  restoring  moment  for 
typical  ships.  For  present  purposes  the  details  of  excitation,  F(t)  need 
not  be  dealt  with  except  for  the  assumptions  that  the  exciting  moment  is 
not  a function  of  roll  angle,  but  is  zero  mean  and  in  the  case  of  random 
excitation  is  Gaussian. 

Thus  the  assumed  rolling  equation  is  written  as: 

..  .2  .3  2 3 

cp+oracp+crcp+(a  /a)<p  + a £cp  = X(t)  (6) 

i 3 

with  X(t)  = F( t )/ I 

and  it  is  assumed,  as  in  previous  work  with  the  same  objectives,  that 
Eq.  (6)  is  a plausible  representation  of  zero  speed  ship  rolling  in  beam 
seas.  Because  of  the  cubic  term  in  roll  velocity  the  previous  objections 
to  the  application  of  the  functional  series  are  removed.  In  their  place 
are  coefficients  of  damping,  ot  , a , which  may  be  chosen  as  functions  of 
the  range  of  variation  of  roll  velocity. 
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THE  FUNCTIONAL  SERIES  MODEL 


The  Time  Domain  Model 


It  will  be  convenient  to  depart  slightly  from  the  previous  notation 
and  consider  the  system  of  concern  to  have  an  "output"  Y(t)  (which  may  be 
cp,  cp,  etc.)  and  an  "input"  X(t)  which  is  assumed  to  be  zero  mean  whether 
it  is  deterministic  or  random.  It  is  assumed  that  the  output,  Y(t)  is  a 
sufficiently  regular  function  so  that  it  may  at  least  be  expanded  in  an 
infinite  functional  series: 


Y(t)=  £ H ’•  Jign(t^ta'-tn)X(t-ti)X(t-t2)-'-X(t-tn)dtidt2...dtnJ  (7) 


(Omission  of  limits  on  integrals  here  and  throughout  this  report  signify 
limits  of  -00  and  •+». ) 


In  Eq.  (7)  the  kernels,  gn(t  •••),  are  "time  invariant"  since  they 


are  assumed  to  be  functions  only  of  time  di fferences.  In  the  present 
application  only  the  input,  X(t)  varies  with  time,  and  the  dynamic  pro- 
perties of  the  system  (the  ship)  are  thus  imbedded  in  the  kernels.  Without 
loss  of  general i ty'^' ^ the  kernels  are  assumed  to  be  real  and  completely 
symmetrical  in  their  arguments.  That  is: 


9n(t  •“O  = 9n(to,t  ,“*tn,t  ) = 
n 1 a n 023  n 1 


for  any  rearrangement  of  the  t^ 


So  long  as  the  series,  Eq.  (7),  converges, i ts  value  is  that  of  the 


kernel  of  zeroth  degree  (gQ)  when  the  input,  X(t)  is  zero.  In  the  present 


application  to  ship  rolling  it  may  be  taken  that  there  will  be  no  roll 
unless  there  is  excitation,  so  that  the  first  term  in  Eq.  (7)  will  be 
dropped  in  all  subsequent  development. 


So  long  as  the  sum  of  the  integrals  of  the  absolute  values  of  all 
kernels  is  finite,  the  series  converges  for  bounded  input,  and  for  stochastic 
input  if  in  addition  the  input  is  strictly  stationary  with  bounded  moments 
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of  all  orders.  These  restrictions  seem  acceptable  physically,  expecially 
since  they  have  been  accepted  for  some  time  in  seakeeping  research  on 
linear  processes. 

It  can  easily  be  imagined  from  the  form  of  Eq.  (7)  that  the  com- 
plexity of  any  answer  which  might  result  would  increase  geometrically  with 
the  order  of  the  term.  Consequently,  it  is  hoped  that  the  non-linearities 
in  any  given  system  are  weak  enough  that  the  series  may  be  truncated  at 
a relatively  few  terms,  in  which  case  it  is  termed  a "functional  polynomial." 

Impulse  and  Frequency  Response  Functions 

The  kernels  in  Eq.  (7)  may  be  considered  as  describing  the  system 
through  a series  of  nth  degree  impulse  response  functions.  It  is  presumed 
that  each  impulse  response  function  is  sufficiently  smooth  and  integrable 
so  that  there  is  no  trouble  about  existence  of  an  n-fold  Fourier  transform. 
Accordingly,  it  is  assumed  that  to  each  nth  degree  impulse  response  func- 
tion there  corresponds  an  nth  degree  frequency  response  function, 

Gn(u>  ,u)^..u>n).  The  transform  pairs  relating  impulse  and  frequency  response 
functions  may  be  defined  as  follows: 

g ( t , t • • • t ) = rf-*f  G (co,ou**aj)  Expi  i E u).t.  dw  dm  •••dui  (6) 

n 1 2 n (2n)n  n 1 2 n L j = i J JJ  1 2 n ' ' 


G ( uu  , <u  • • • tu  ) = I • f g (t  ,t  • • • t ) Exp!  - i £ U).  t . I dt  dt  • • • dt 

n 1 2 n “ “ " n'  1 2 n'  rL  J JJ  1 2 n 

The  first  degree  frequency  response  function  is  the  familiar  linear 
one.  The  second  degree  frequency  response  function  is  the  one  treated  in 
Refs.  15,  16,  and  17  in  an  application  to  added  resistance.  Regardless 
of  the  degree,  the  basic  importance  of  the  transform  of  the  impulse  re- 
sponse function  is  the  same;  that  is,  convolution  in  the  time  domain 
corresponds  to  multiplication  in  the  frequency  domain.  In  practical 
applications  it  is  most  convenient  to  work  in  the  frequency  domain. 

As  a consequence  of  the  assumed  symmetry  of  the  impulse  response 
functions  and  the  transform,  Eq.  (9),  the  nth  degree  frequency  response 
function  is  also  symmetric  in  its  arguments.  That  is, 
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G (u>,U). ..<!))  = G (u)  ,0)  . ..(!))  ... 
n'  x a n7  n'  27  1 n7 


for  any  and  all  rearrangements  of  the  u>  . Additionally,  because  the 
inpulse  response  functions  are  real: 

G (-U), -co, •••-(«)  = G (00,(0.. .001 
n'  17  37  n7  nv  x7  2 n7 

where  the  star  denotes  the  complex  conjugate,  and  a 1 1 arguments  on  the 
left  hand  side  are  negative. 

Linear  Operations  on  Input  and  Output 

A knowledge  of  the  effects  on  the  various  frequency  response 

functions  of  linear  operations  on  input  and/or  output  is  useful  for  both 

analyses  and  interpretation.  Figure  2 indicates  the  type  of  cascading 

assumed.  (Theory  for  more  complicated  cases  may  be  found  in  Barrett'0 
20* 

and  George  .)  In  the  middle  of  the  block  diagram,  the  series,  Eq.  (7), 
is  indicated.  Each  term  of  the  series  is  excited  by  X(t)  and  the  result 
for  the  nth  term  is  noted  as  Yn( t ) . The  Yn(t)  are  summed  to  produce  Y(t). 
A linear  input  filter  acts  upon  W(t)  to  produce  X(t),  and  this  filter  is 
characterized  by  impulse  response  function,  4(t).  Similarly,  at  the  other 
end  of  things,  the  output,  Y(t)  is  acted  upon  by  another  linear  filter 
(with  impulse  response  function  h(T))  to  produce  a new  output,  Z(t). 

Because  both  input  and  output  filters  are  linear  it  is  necessary 
only  to  consider  the  influence  of  filtering  on  the  nth  term  of  the  series: 


Yn(t)  = 9n(t17t2"tn)  x(  t-O*  *X(  t-t^dt^  • . dtp  (10) 

Considering  the  input  filter  first: 

X(t-tj)  = J *(t.)  W(t-tj-T.)dT.  (U) 

Then  substituting  Eq.  (11)  in  Eq.  ( 10),  making  a change  in  variable  of 
the  form  v.  = t^  + t and  rearranging: 


20.  George,  D.A.,  "Continuous  Non-Linear  Systems,"  Doctoral  Dissertation, 
Department  of  Electrical  Engineering,  M. I.T.,  July  1 9 59 - 
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where: 


v ...v  ) W(t-v  ) W(t-v  )***W(t-v  )dv  dv  '-dv 
2 rr  j v 27  ' n'  a 2 n 


(12) 


q^i  v , v . • • • v ) = rr**  r g ( v - t ,v-t  ,*  • v - t ) 

sn  v S 2’  n7  JJ  J yn'  1 ^ 2 2 n n7 


£(t  ) 1(t  )•••£(  t )dT  ,dT  ...dT  (13) 

1 2 n 1 2 n 

Noting  that  Eq.  (12)  is  of  the  same  form  as  Eq.  ( 10),  Eq.  (13)  is 

an  expression  for  the  n degree  impulse  response  function  relating  input 

W(t)  to  output  Y (t).  Substituting  for  g (v  -T  ,*••)  in  the  right  hand 
n n 1 1 

side  of  Eq.  (13)  its  transform,  Eq.  (8),  and  noting  that  the  filter  fre- 
quency response  function  is: 

. J dT 

1-  n -| 

• G ( 12  , uu  • . a)  ) Expi  i £ a),  v.  dv  • • • dv  (lU) 

n'  i7  2 n7  KL  j = , j jj  1 n v 7 

From  which  the  nth  degree  frequency  response  function  relating 
input  W(t)  and  output  Y ( t)  is  found  from  Eq.  (8): 


L(») 

there  results: 


WY,  \ 

g ( V , V , . . . V ) 
an  v i 7 27  n7 


WY 

Gn  (“  »“a»  ••*%)»  L(®)L( «i>a)**L(o>n)Gn(o)  ,«Ba*-oun)  (15) 

Now  considering  the  output  filter,  the  component  of  output  corres- 
ponding to  the  nth  term  of  the  series  is: 

Zn(t)  = J h(T)  Yn(t-T)dT  (16) 


Substituting  Eq.  ( 10)  in  Eq.  ( 16),  and  after  making  a variable  change  of 

the  form  t = t + t,  there  results: 

1 l 
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2n<0  “ JH  S^VY'V  X(t‘Tx)  X(t-Ta)...X(t-Tn)dT  ...dTn  (17) 

where 

9JZ(Tl,Y*Tn)  = I h(T)  9n(Tx-T>YT,‘"VT)dT  O8) 

Again  noting  the  similarity  of  form  between  Eq.  (17)  and  Eq.  (10),  Eq.  ( 18) 
is  an  expression  for  the  nth  degree  impulse  response  function  relating 
input  X(t)  to  filtered  output  Zn(t). 

As  before,  substituting  for  g (t  -t  •••)  in  the  right  hand  side  of 

n 1 

Eq.  (18),  its  transform,  Eq.  (8),  and  noting  that  the  output  filter  fre- 
quency response  function  is: 

H(  id)  = f h(  T)e"  1 ^ dT 

There  results  an  expression  of  the  form  of  Eq.  (8)  from  which  the  nth 
degree  frequency  response  function  relating  input  X(t)  with  filtered 
output  Zn(t)  is  found  to  be: 


XZ  / n \ 

3 (id,u>,**id)  = G(u;,id««id)H\E  u >.  ) 

n 1 2'  n n'  x 2 n'  j / 


Combining  Eq.  (15)  and  ( I9)  with  reference  to  Figure  2,  the  nth 
degree  frequency  response  function  relating  input  W(t)  to  filtered  output 
component  Zn(t)  becomes: 


„WZ,  % „WY 

G ( (D  , CD  • • (D  ) = G 
n ' 1 2 n ' n 


GWY(  id  , u>  • • • uo  ) h ^ £ id.  ) 
n i 2 n \ j * 

/ n 

L(«  )L(oj  )--L(c»  )G  (»,«••«>  ) Hi  S u). 

1 2 n n 1 2 n ' j_j  j 


7 


R- 1 894 


DETERMINATION  OF  FREQUENCY  RESPONSE  FUNCTIONS 
CORRESPONDING  TO  THE  ROLL  EQUATION 


The  general  objective  of  the  next  step  in  application  of  the  func- 
tional series  is  to  expand  the  equation  in  the  series  and  to  develop 
expressions  for  the  frequency  response  functions  in  terms  of  coefficients 
in  the  equation.  Because  a series  is  involved  at  the  outset,  a series  of 
frequency  response  functions  can  be  expected  to  result.  To  gain  some 
appreciation  as  to  how  high  a degree  the  series  might  be  carried,  it  is 
instructive  to  consider  the  result  when  the  series,  Eq.  (7)  is  excited 
by  a simple  harmonic  function. 

The  excitation  is  assumed  in  the  form: 

X(t)  = X cosiut 
' ' o 


X 1 r 

= ~ £ Exp[_(-l)r  iu)t 


. S r 


Substitution  of  this  excitation  in  the  general  n*  term  of  Eq.  (7);  and 
after  some  manipulation  and  the  application  of  Eq.  (9),  the  result  is: 
x"  1 1 1 , r,  r2  r v 

Yn(t)="H  2 2 •**  2 Gn  l(-0  <**(-U  «V(-0  V* 

2 ri=o  r2=o  rn=o 

r~  n r . -j 

ExpLi  aJt  Z(-l)Jj  (22) 

j = l 

When  ships  or  models  are  experimentally  subjected  to  beam  regular 
waves  of  frequency  UJ,  the  dominant  part  of  the  roll  response  is  at  this 
frequency.  Higher  harmonics  or  shifts  in  the  mean,  if  detectable,  are  of 
much  smaller  magnitude  than  the  response  at  frequency  uu.  When  present, 
the  nonlinearity  manifests  itself  as  a departure  from  linearity  with  wave 
height  of  the  apparent  rolling  amplitude.  Accordingly,  the  dominant 
frequency  response  functions  would  be  expected  to  be  those  which  make  a 
contribution  to  the  amplitude  of  the  component  of  response  at  frequency  uu. 


R_ 1894 


Considering  the  expression  for  the  contibution  to  the  n degree 
term  to  the  output,  a contribution  to  response  at  frequency  uj  happens  only 


■ 

I 


E (-1)  j = ±1 

j=« 

in  Eq.  (22). 

Noting  that  r^.  is  either  zero  or  one  it  may  be  seen  that  no  combina- 
tions of  the  n in  Eq.  (22)  can  produce  the  above  result  if  n is  even.  Thus 
none  of  the  terms  of  even  degree  in  the  series,  Eq.  (j)  contribute  to 
response  at  frequency  u),  and  accordingly  it  must  be  anticipated  that  none 
will  be  dominant  in  the  present  application. 

When  n is  odd  in  Eq.  (22)  a contribution  to  output  at  frequency  u; 
occurs  when  (n+l)/2  of  the  r.  are  zero.  There  are 


/ n+K  ' / n- 1 \ ' 

( -5— ; • I-"—  J • 

such  terms,  and  for  each  the  argument  of  the  exponential  is  ( i cut) . Because 
of  the  symmetry  of  the  frequency  response  functions  all  of  these  terms 
are  equal.  In  addition,  there  are  the  same  number  of  terms  (all  equal  by 
symmetry)  for  which  the  argument  of  the  exponential  is  ( — i out)  and  in  which 
(n-l)/2  of  the  r.  are  zero,  so  that  this  latter  contribution  is  the  complex 
conjugate  of  the  first. 

Accordingly,  the  total  contribution  of  the  nth  degree  term  in  Eq.  (7) 
to  the  response  at  frequency  u)  becomes  (for  n odd) : 


VR  n ' 

X n • 
o 


/ n+1  x » / n *“  1 \ t 


Re  [ Gn(  &,  uu,  - a),  uj,  - • • ) J cos  (JJt 

- lm[Gn(  uu,  a), -uj,  (u,-ou,  • • ) ] s inoit 


Evaluating  for  odd  n,  the  response  at  frequency  u)  of  the  system 
defined  by  Eq.  (7)  to  excitation  of  the  form  X(t)  = Xq  cosuut  is  itself 
of  the  form: 
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P cosiut  - Q.  sinu)t 

where: 

P = X RetG  (uj)] 

0 1 

+ ^ X^  Re  [G3(  u^-uj)  j 
+ H Xo  R6^G  ( u),id,-u),id,-u))] 

Or 

+ Xo  Re^7(tu,  u),-tu,u),- cm)  ] 

+ •••  (24) 

and  Q is  the  same  except  that  the  imaginary  parts  of  the  frequency  responses 
are  involved. 

It  is  clear  from  this  result  that  the  series  must  be  carried  at 
least  to  third  degree  in  order  to  reflect  the  expected  type  of  nonlinearity, 
and  if  this  is  not  sufficient,  then  to  fifth  degree.  For  present  purposes 
it  was  assumed  that  carrying  derivations  through  to  fifth  degree  would  be 
more  than  sufficient. 

To  proceed  it  is  convenient  to  consider  an  equation  equivalent  in 
form  to  Eq.  (6),  the  roll  equation: 

A Y(t)  + A (Y(t))3  + B Y(t)  + B (Y(t))3 

13  13 

+ C Y(t)  + C (Y(t))3  = X(t)  (25) 

1 3 

In  the  above, Y(t)  corresponds  to  roll,  or  output;  X(t)  is  excitation, 
and  it  is  assumed  that  the  equation  is  stable,  possesses  one  solution, 
and  that  Y(t)  vanishes  when  X(t)  does.  It  is  presumed  that  the  coeffi- 
cients with  subscript  1 are  non-zero,  but  that  any  or  all  of  the  coeffi- 
cients with  subscript  3 may  be  zero. 

The  approach  to  the  expansion  is  that  called  the  "harmonic  input 
method"  by  Bedrosian  and  Rice.''  In  this  method  the  input,  X(t),  is 
assumed  to  be  the  summation  of  a very  large  number  of  harmonic  functions 
with  incommensurate  frequencies,  say: 
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X(t)  = £ Exp(it(u.)  (26) 

J=1 

Substituting  Eq.  (26)  into  Eq.  (7),  omitting  the  constant  term,  and  after 
the  application  of  Eq.  (9),  there  results: 

Y(t)  = I G (uij)  Exp(  icujt) 

J=1 

+ E EG  (u)  ExpLit(  u>n  + u>  .0)] 

Jl=l  J2=l  2 J1  J2  J 

+ E Z ^ G3(u;jl’UJj2  ,a;J3^  Exp^‘t^U)Jl  + UJJ2+U3J3  ^ 

J 1 — 1 J2= 1 J 3~ • 

+ Jl-l  "jf-l  G4(U)j1>‘’'UJj4)  ExpCit(U)Jl+aj2+UJJ3+U)J4)] 

+ E •••  E G (ujiu‘,'U)is)  ExpCit(u)n+**  + u)  c)^l 
Jl=l  J5=l  5 J 


For  convenience  two  notations  may  be  established: 

"•  "jn>  * HnO’V--") 

Mi 1(0,,+  •••  «'Jn)j  - En 

Now  considering  the  last  two  terms  on  the  left  of  Eq.  (25),  and  substituting 
Eq.  (27)  and  renumbering  indices  where  required: 

C Y(t)  + C (Y( t))3  = I C H (1)E 
1 3 Jl  = l 1 1 1 

+ E E C H ( 1 , 2 ) E 
J 1 = 1 J2= 1 1 2 2 

+ E •••  E [C  H (1,2,3)  + C H ( 1 )H  (2)H  ( 3) 3 E 

J 1= 1 J3=l  13  3111  3 

+ E • • • E [C  H (1,2,3,^)  + 3C  H (1)H  (2)H  (3,'OjE 
Jl=l  jU=|  1 4 3113 


I > 


Ejjrvii 


■ E [C  H (1,2,3A5)+3C  H (1)H  (2,3)H  (U,5) 
J 1=1  J5=l  1 6 312  2 


+ E 


+ 3C  H ( 1 )H  (2)h  (3^,5)]  E 

3 113  b 


+ • 


(28) 


Considering  the  differentiation  of  Eq.  (27)  with  respect  to  time, 
the  general  term  in  Y(t)  becomes: 

E • 


E i ( u)  . , + ••  + (!),  ) G ( u)  . , , * * • . ) E 

..  , , . ' Jl  Jn'  ny  jr  Jn  n 

Jl  = l Jn= 1 


and  similarly  the  typical  term  in  Y(t)  is: 


E •••  E -(  u)  +•  . + 0) . )2  G (tu. )E 
, J1  Jn7  n'  j 1 ’ Jn'  n 

J 1 = 1 Jn=  1 


Accordingly,  the  expansion  of  the  velocity  terms  in  Eq.  (25)  becomes 


B Y(t)+B  (Y(  t))3  = E iui  B H ( 1)E 

J 1 1 1 1 


Jl  = l 


+ E E i(u»  +u»  )B  H (1,2)E 
Jl-1  J2=  1 J1  J2  2 2 


+ • • • 


(29) 


and  that  of  the  acceleration  terms: 


A Y(t)+A  (Y(t))3  = E -ttj.A  H (l)E 
1 3 ji=i  Jl  1 1 1 


+ E 
Jl  = l 


E 

J2=  1 


" (U)ii+u>i?)2a  H ( '»2)E 
J 1 ocL  12  2 


+ • • • 


(30) 


Completing  the  expansions  Eq.  (29)  and  ( 30 ) through  terms  in  E and  adding 
the  results  to  Eq.  (28),  the  expansion  of  the  left  hand  side  of  Eq.  (25) 
through  terms  of  fifth  degree  is  completed.  The  result  may  be  written 
as  follows: 
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X(t)  = £ E = (left  hand  side  of  Eq.  (25)) 

ji=i  1 

= Z Q E 
Jl  = l 1 * 

+ Z Z Q E 
Jl=l  J2=l  2 2 

+ Z • • • Z Q E 
Jl=l  J3= 1 3 3 

+ Z ...  I Q E 
Jl  = l j4=l  4 4 

+ Z •'*  z Q E 
Jl  = l J5=l  5 5 


Where  the  Q's  are  defined  as  follows: 

\ - Dx < 1 ^ji > vo 

Q2  - D1<,<®J1+0,J2»  Hs(1’2) 

Q3  = Di(l(a)J1+ ««J2+ a,J3))  H3(  1,2,3) 

+ D3(-iu)jrj2a’j3)  \{])  Hi(2)  Hx(3) 

Q4  = D (l(«ttJ1  + «,J2+  -*^)  H4(1,2,3,M 

+ 303(-i^j,  ®j2("j3+aju))  ^<0  h2(3,4) 

Q&  = Di(1(u;j1+u>j2+..+uJj5))  Ml, 2,3, 4, 5) 

+ 3D3(-iu,Jl(u,J2+U,J3)(u,j4+u,J5>)  H1(0  Hs(2'3)  H2(^5) 
+ 3D3^"iuJl  aJ2^  U>J3+<UJ4+U,J5))  « (0  V2)  H (3,4,5) 


and  where: 


D (0)  = A £V  + B 0 + C 
n n n n 


Because  the  frequencies  in  the  excitation,  the  uk,  are  assumed  to 


be  incommensurate,  Eq.  (31)  can  be  satisfied  through  terms  of  fifth  degree 
only  if: 


E E = £ Q E 

jl=l  1 Jl=l  1 1 


0 = I E Q E 


Jl=l  J2=l 


0 = £ • • £ H E 

jl=l  J3=l  3 3 


0 = £ £ HE 

Jl=l  jM  4 4 


0 = £ •••  £ Q.  E 

Jl=1  J 5=  1 6 5 


Essentially  what  is  assumed  is  that  none  of  the  possible  exponentials  in 
En  is  the  same  as  any  of  the  exponentials  in  En+j,  etc.  Each  equation  of 
degree  (n)  is  satisfied  by  equating  to  zero  or  one  all  the  terms  in 
which  have  the  same  exponential  time  factor  for  an  arbitrary  choice  of  n 
of  the  frequencies  in  the  input. 


Considering  Eq.  (38)  and  noting  that  E = Exp(itu)jj),  it  is  neces- 


sary that  Q.  = 1 for  any  choice  of  J1  (and  thus  any  particular  frequency), 
l 

Thus  from  Eq.  (32): 


D ( i a)  G ( u>)  = 1 
1 x 


G ( uu)  = 

i D ( i oj) 

l 


= 1/(_oFa  +i<JJB  +C  ) 

l l l 


This  is  the  result  expected  if  all  the  coefficients  in  the  cubic  terms  of 
Eq.  (25)  are  zero. 


Now  considering  Eq.  (39)  and  making  the  substitution  for  Q^: 


°-  j£,  D,<,<<VV)  G2<  W’  ExpOtfw 
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Any  two  of  the  incommensurate  frequencies,  say  <d  and  id  , result  in  a unique 

1 a 

time  factor,  it(u)+<D  ),  and  there  will  be  two  terms  in  the  double  sum  which 
1 2 

involve  this  same  time  factor.  Thus  the  equation  is  satisfied  if  for  every 
pair  of  frequencies: 

0 = D ( i(  u)  + <d  ))  G ( <d  , (D  ) + D ( i(  id  + id  ))  G ( id  ,<d  ) 

1 13  213  1 31  331 

Since  0 (•••)  is  not  zero  and  the  frequency  response  functions  are  symmetric: 

1 


G ( <u  ,u>  ) = 0 

3 13 


m 


Equation  (Uo)  for  the  third  degree  frequency  response  function  may 
be  written: 


0 = £ 


T ; 0 (l(tt>ll+«10+«J)  G (<•»,„»  o,«ia) 


Jl  = l J2= 1 J3=l 


Jl  J2  J3"  V Jl’  J2  J3' 


* W VV  G1<“’J3>J  ' 

Again,  any  three  of  the  incommensurate  frequencies  result  in  a unique  time 
factor.  For  any  choice  of  three  frequencies  there  are  six  terms  in  the 
triple  sum  with  the  same  time  factor,  and  these  terms  correspond  to  all 
the  permutations  in  order  of  the  three  frequencies.  Because  all  terms 
are  symmetric  in  the  equation,  all  six  terms  are  equal. 

Thus  the  equation  is  satisfied  if  for  any  particular  choice  of  three 
frequencies,  say  id  , id  , <d  ; 

12  3 

0 = D ( i(  w +w  +w  ))  G ( w , w , <i)  ) 

1 133  3123 

+ 0 (-iuucuuu)  G (u))  G (00)  G (id) 

3 123  11  13  13 


and  thus: 


G (id,id,(d)  = -D  ( - i id  u>  (D  ) G (id+id+uu) 

3123  3 123  1123 


G(u>)G(<d)g(<d) 

11  12  13 


(45) 
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The  relationship  to  the  coefficients  in  the  differential  equation  is  pro- 
vided by  D (•••)  and  the  previously  obtained  relationship  for  G ( u,),  Eq. 

3 l 

(43). 

Taking  advantage  of  Eq.  (44)  for  the  quadratic  frequency  response, 
Eq.  (41 ) may  be  written: 

°-  V'WW  • 

J 1 = I J 4=  1 

Cxp[tt(»JI  + «B+  ••  *Jt)J 

The  solution  in  this  case  is  essentially  the  same  as  that  for  G (cu.ou  ). 

2 12 

For  every  choice  of  four  frequencies  which  define  a particular  time  factor 

there  are  24  terms  which  have  the  same  time  factor.  These  correspond  to 

all  permutations  in  order  of  the  frequencies,  and  because  of  the  symmetry 

of  all  factors  of  each  term,  all  2k  terms  are  equal.  Again  since  D (•••) 

1 

is  not  zero  there  results: 

G ( u>  , (u  ,uj  , (u  ) = 0 (46) 

4 12  3 4 

The  solution  for  the  fifth  degree  frequency  response  comes  from 

Eq.  (42).  Taking  advantage  of  the  relationship  for  G ( u>  , 00 ) in  Eq.  (44), 

212 

Eq.  (42)  becomes: 


+ 3D3<-i“jlV“j3+"jt+'"j5))  W 

• G^iDjsj)  Exp[it(u)J1+.  .+»J5>] 

In  the  above  for  a choice  of  five  of  the  incommensurate  frequencies  which 
result  in  a particular  value  of  the  time  factor,  there  are  120  terms  in 
the  coefficient  of  the  exponential  which  have  the  same  form  as  the  general 
term  above.  These  correspond  to  all  the  possible  permutations  in  frequency 
order.  Both  factors  of  the  first  part  of  the  general  term  above  are 
symmetrical  so  that  its  contribution  to  the  total  coefficient  of  the 
exponential  is  straightforward.  The  second  part  of  the  general  term  is 
not  completely  symmetrical,  its  contribution  works  out  to  10  groups  of 
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12  identical  terms  each.  As  in  the  previous  developments  the  coefficient 

of  the  exponential  for  an  arbitrary  choice  of  5 frequencies  ( u>  , • • • id  ) 

was  equated  to  zero  and  this  equation  solved  for  G (<d  •••).  Since  the 

5 1 

result  was  a function  of  G ( id  • • • ) as  well  as  G (u>),the  third  degree 

3 1 1 

frequency  response  function  was  eliminated  by  use  of  Eq.  (45)  and  the 
final  result  is  given  in  Eq.  (47). 

G ( id  , u>  • • id  ) = ~G  (uj+id+id+id+id)  * 

6 1 2 5 10  1 1 2 3 4 6 


where: 


G (id  ) G (id)  G (id  ) G (id)  G (id)  K (<d  , u>  , id  , id  , id  ) (hj) 

11  IS  13  14  15  512345 


K(u>,u>,id,«d,u>)  = K(  12345)  + K(  13245)  + K(  14235) 

5 1 2 3 4 5 

+ K( 15234)  + K(23l45)  + K(24l35) 

+ K( 2 5 1 34 ) + K( 34125)  + K( 35124) 

+ K(45123)  (48a) 

and: 


K(jkmnp)  = °3  ( " ' “j- 


U)  +ID  +(JD  )) 
m n p" 


0 ( - i U)  id  id  ) 
3 m n p 


G ( id  +u>  +(D  ) 
1 m n p 


(48b) 


To  summarize  the  section,  the  first  five  frequency  response  functions 
resulting  from  the  expansion  of  Eq.  (25)  are  given  in  Eqs.  (43)  through 
(48).  All  the  functions  corresponding  to  even  degrees  are  zero,  a result 
consistent  with  the  results  of  the  initial  development  in  this  section. 

The  first  degree  function  is  just  what  would  have  been  obtained  had  there 
been  no  cubic  terms  in  Eq.  (25).  The  third  degree  function  looks  not  too 
complicated  but  the  fifth  degree  function  is  quite  cumbersome. 
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SCALAR  SPECTRUM  OF  RESPONSE 
OF  THE  FUNCTIONAL  SERIES  MODEL 

The  main  concern  herein  is  with  random  excitation.  For  this  purpose 
it  will  be  assumed  that  the  excitation,  X(t),  in  Eq.  (7)  is  a stationary, 
Gaussian,  zero  mean  process.  The  autocorrelation  of  the  process  may  be 
denoted: 


Rx(  t)  = X(t)  X(t-T) 

The  two-sided  spectrum  of  the  process  is  by  definition: 
Sx(u>)  = J R (t)  ExpC-iWTjdT 


The  integral  over  positive  and  negative  frequencies  of  the  spectrum, 
(Eq.50),  is  2nRx(0);  that  is,  2tt  times  the  process  variance. 


Bedrosian  and  Rice  develop  a series  for  the  two-sided  scalar 
spectrum  of  response,  Sy(uo),  given  such  excitation  and  the  functional 
series  model,  Eq.  (7),  with  the  constant  term  omitted.  Their  result  was 
expanded  to  include  all  contributions  of  terms  in  the  series  up  to  the 
fifth  degree;  that  is,  the  expansion  was  carried  out  as  though  the  series, 
Eq.  (7)  had  been  truncated  after  the  term  involving  the  fifth  degree 
impulse  response.  The  general  result  for  the  scaiar  spectrum  of  response 
(in  present  notation)  is  given  in  Eq.  (51),  with  further  definitions  in 
Eq.  (52)  through  (56). 

Sy(“)  = SY|(  UJ)  + Sy2(^  + Sy3(u>)  + + SY5^^  + " (50 

where  the  contributions  to  Sy(  uu)  are  as  follows: 
syiH  = SXH  ^ I Ga(u),  v,-v)  Sx(v)dv 

2 

+ 5 JJ  G ( u),v,-v,w,-w)  S (v)  S (w)dvdw  + ••  (52) 

(2v)  5 xx 
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SY2^^  = 7t  I sx(““  u)  Sx(u)  Gg(®-M»u) 

2 

+ 4 I g4(U)'u>u> v>_v)  sx(v)dv  + “ du  (53) 

sy3(uj)  = ~ J.f  SX((U“U-V)  sx(u)  sx(v)  G3(uj-u-v,u,v) 

2 

+ J G ( uu-u-v,  u,v,w,-w)  S^(w)  dw  + • • dudv  (5*+) 

SY^(  JJJ  sx(u>“u”V-w)  Sx(u)  Sx(v)  Sx(w)  • 

2 

G ( u>-u-v-w,  u,  v, w)  + ••  dudvdw  (55) 

4 

Sy^(  “0  = ~~  J’J’JJ’  Sx(w-s-v-u-w)  Sx(s)  Sx(v)  Sx(u)  Sx(w)  • 

2 

G ( u)-s-v-u-w,s,v,u,w)  + **  dsdvdudw  (56) 

6 

As  can  be  noted,  each  term  in  the  series  is  itself  a series.  The 
series  for  Eqs.  (52),  (5*0  and  (56)  involve  only  the  frequency  response 
functions  of  odd  degree,  and  the  series  for  Eqs.  (53)  and  (55)  involve 
only  the  functions  of  even  degree.  The  first  and  second  degree  frequency 
response  functions  appear  explicitly  only  once  in  Eq.  (52)  through  (56), 
functions  of  the  third  and  fourth  degree  only  twice,  and  the  fifth  degree 
function  appears  three  times.  The  form  of  the  series  is  such  that  con- 
tinuation of  the  expansion  to  terms  of  higher  degree  involves  no  new 
occurrences  of  the  first  through  fifth  degree  frequency  response  functions. 

The  foregoing  is  the  generalized  scalar  spectrum  of  output,  Y(t), 
of  Eq.  (?)•  The  relationships  between  output  spectrum  and  the  spectra  of 
output  velocity  and  acceleration  are  also  of  interest  in  the  present  ap- 
plication. The  differentiation  process  may  be  considered  a linear  operator 
operating  upon  Y(t).  Considering  this  operator  as  an  output  filter. 

Figure  2,  the  nth  degree  frequency  response  function  relating  output 
velocity,  Y(t)  to  input,  X(t)  becomes  (from  Eq.  19,  replacing  H(o')  by  io): 
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G*^  ( 0)  , • • • u)  ) = i GX^(  u>  • • • ) £ u) 


and  that  relating  output  acceleration  to  input  is: 

2 

GnV  < = " G^y  (Ui‘‘*)  (A  “j)  (58) 

Because  the  relations  Eq.  (50  through  (56)  are  for  generalized 
output,  the  spectrum  of  velocity  would  result  if  each  frequency  response 
function  in  Eq.  (51)  through  (56)  were  replaced  by  that  relating  ouput 
velocity  to  input,  Eq.  (57)-  It  happens  that  the  sum  of  the  frequency 
arguments  of  the  portion  of  every  frequency  response  function  required  in 
tq.  (52)  through  (56)  is  equal  to  oj.  Accordingly,  the  expression  for 
velocity  spectrum  differs  from  that  of  Eq.  (51)  through  (56)  only  by  a 
factor  of  ( i ou)  inside  the  squared  absolute  values,  and  thus  the  spectrum 
of  output  velocity  becomes: 

S f(  <«)  = Sy(  u>)  (59] 


The  same  argument  applies  to  accelerations  so  that  the  spectrum  of  output 
acceleration  becomes: 

Sy(  <“)  = w4  Sy(  u>)  (60) 

Thus  the  relationships  between  output,  and  output  acceleration  and 
velocity  spectra  for  the  non-linear  model  are  the  same  as  those  used  in 
the  1 inear  context. 
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SCALAR  SPECTRUM  OF  SINGLE 
DEGREE-OF-FREEOOM  ROLL  RESPONSE 

To  evaluate  the  scalar  spectrum  of  single  degree-of-f reedom  roll 
response  the  expressions  for  the  frequency  response  functions,  Eq.  (1+3) 
through  (1+8),  must  be  substituted  into  Eq.  (50  through  (56).  The  result 
of  this  operation  would  be  the  spectrum  of  response  of  the  system  defined 
by  Eq.  (25)  in  terms  of  the  coefficients  in  the  equation.  Conversion  to 
the  spectrum  of  roll  response  would  then  be  achieved  by  replacing  these 
coefficients  with  the  corresponding  coefficients  in  Eq.  (6)  which  is  the 
assumed  rolling  equation. 

It  is  useful  for  later  purposes  to  consider  the  result  when  all 

the  coefficients  of  the  cubic  terms  in  Eq.  (6)  or  (25)  are  zero,  that  is: 

the  equations  are  totally  linear.  In  this  case  it  can  be  seen  from  Eq.  (43) 

through  (48)  that  the  only  non-zero  frequency  response  function  would  be 

G (a).  Thus  the  general  expression  for  the  spectrum  (Eq.  51  to  56)  breaks 
l 

down  to  the  usual  relation  between  input  and  output  spectra  for  linear 
systems . 

The  "linear  spectrum"  of  roll  output  (the  spectrum  if  the  coeffi- 
cients of  nonlinear  terms  are  zero)  may  be  defined  as: 

- |g  («o|a  s («)  (6i) 

T 1 1 I ^ 


also,  the  "linear"  spectrum  of  roll  output  velocity  may  be  defined: 


Y.(  (t)  = (JU 
cp'  ' 


G (u>) 


S (a) 


(62) 


Both  these  spectra  are  two-sided,  as  is  S (it). 


Continuing  with  definitions,  the  variance  of  1 i near  output  (V^)  and  the 
variance  of  I i near  output  velocity  (V^)  will  be  defined  to  be: 


\ - is  J V “)d 


U! 


(63) 
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V9  " 5n  J Yi  < “>d®  <6h' 

Now  substituting  the  expressions  for  the  nonlinear  frequency  re- 
sponse functions  (Eq.  44  through  47)  in  the  expression  for  the  general 
output  spectrum,  Eq.  (50  through  (56),  and  taking  account  of  the  defini- 
tion, Eq.  (6l ),  the  spectrum  of  nonlinear  roll  response  expanded  through 
terms  in  the  5th  degree  frequency  response  function  may  be  written: 


SqH  = sq,i(°»)  + V «)  + S^w)  + 


where: 


V”9  - y»>  ' I °3(i,“#2>  Vv)dv 

2 

+ G ( id)  JJ  K (u),v,-v,u,-u)  Y (v)  Y (u)dvdu  + •• 

gTT  I B 9 9 

V9  lG1(”)l  II  V"'u"v)  Vu)  Vv)dudv  ■ 

- D ( iuv(u+v- 0)  +—  J K ( ou-u-v,  u,v,  w,-w)  Y^(w)dw  + •• 

v9  -~l  KH  mi  Vs'  Vv)  Vu)  vw)  ■ 


K ( u>-s-v-u-w,s,  u,v,w)+ 


dsdvdudw 


In  Eq.  (66)  through  (68),  G ( id)  is  defined  by  Eq.  (43),  D ( i “>•• ) is  defined 
by  Eq.  (37),  and  K ( <d, ^)  is  defined  in  Eq.  (48). 

5 

Perhaps  the  most  notable  thing  about  the  roll  response  spectrum  as 
shown  by  Eq.  (65)  through  (68)  is  that  the  input  spectrum,  S ( id),  is  no 
longer  an  explicit  part  of  any  of  the  expressions.  It  is  embedded  in  the 
linear  roll  response  spectrum, Eq.  (6l).  In  this  respect  the  present 
development  begins  to  resemble  that  of  Yamanouchi^. 
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In  practice,  estimates  are  required  for  the  realizable  single  sided 
roll  spectrum.  In  order  to  convert  the  expressions  in  Eq.  (66)  through 
(68),  each  doublesided  spectrum  must  be  replaced  by  ( n)  times  the  corres- 
ponding single  sided  spectrum  with  an  absolute  value  in  the  frequency 
argument.  This  results  in  spectrum  area  equal  to  variance.  To  be  specific 
the  linear  roll  spectrum,  t ( id)  will  be  replaced  by  tt  u ( I u>  I ) where  U(  <d) 
is  the  single  sided  linear  rolling  spectrum  defined  for  positive  frequency. 
Similarly,  each  component  on  the  left  side  of  Eq.  (66)  through  (68)  may  be 
modified  in  the  same  way;  that  is,  S^j(uj)  will  be  replaced  by  ™Ur  j ( l a | ) , 


Comparing  Eq.  (6),  (25)  and  (37)>  the  function  D ( iH)  may  be  written 

3 

D ( ift)  = + ifx*  /a  (69) 

3 3' 

Also,  from  Eq.  (1+3)  and  (6): 

G ( uc)  = l/Lo^-uF  + i oja  aj  (70) 

1 1 

a 2 

!g  (id)I  = l/Uo2-^)2  + ( ux*  0)  j (71) 

1 1 

Now  applying  these  relations  as  far  as  is  possible  analytically, 
the  single  sided  roll  spectrum  (corresponding  to  the  representation  of 
roll  by  Eq.  (6))  may  be  written 


UqH  = 11^(0))  + U93(uj)  + U95(<d)  + ••• 


where: 


U . ( uu)  = U(  (Jo)  G ( id)  Ccr2(  l-3CV,„)-uf  ] + i uuo  o(  l-3a  v./o  cr2) 

H'l  1 x 1 3 H7  i 


+ ^ J J K ( <ju,v,-v,u,-u)  U(v)u(u)c’vdu  +• 


2 " « 5 

o o 


2 


U^(  <o)  = J |G  («)|  JJ1  dvdu  U(  I OJ-u— v I ) U(lul)  U(lvl)  • 

00 

- Co2  + i»  uv(u)-v-u)/c?  + 3 J K ( to-u- v,  u,  v,w,-w)  U(w)dw  +•• 
3 o 5 

2 

U^Coj)  = 15  |G  (u>)l  JJ/J  U(  I od-s-v-u-wI  ) U(lsl)  U(lvl)  U(lul)  U(  lw| ) 


m 


K ( CO-S-V-  u-w.s,  u.v.w)  + 

10  s'  * ’ ’ ' 


dsdvdudw 


(75) 


And  the  function  K (••••)  may  be  evaluated  in  terms  of  coefficients  ‘n 


Eq.  (6)  by  means  of  Eq.  (48)  and  the  relations  Eq.  (69)  and  (70). 


In  the  foregoing  expressions  it  may  be  noted  that  only  one  of  the 
integrals  in  Eq.  (66)  through  (68)  could  be  evaluated.  The  integrals 
involving  K (•••)  have  resisted  all  attempts  at  simplification. 

If  it  is  permissible  to  neglect  the  terms  originating  with  the 
fifth  degree  frequency  response  function,  the  roll  spectrum  becomes: 

10)  = U(  10)^(113)1  | [o2(  l-SCV^-uf]  + (iooa)2  [l-3a  o2]  j 

+ | 1g  (®)l  JJ  U(  l<o-u-vl)  U(  lul)  U(  Ivl  )dudv 

2 

+ | Ig  ( a)  I (a  /of  U(ltu-u-vl)  Cl(lui)  U(lvl)dudv  (76) 

1 3 


where  U(<o)  = the  single  sided  linear  roll  veloc i ty  spectrum  corresponding 
to  Eq.  (62). 


A TRIAL  EVALUATION  OF  EQ.  j6 


In  the  work  of  Ref.  6 an  equation  involving  quadratic  damping  was 
simulated  in  the  time  domain  and  spectra  were  estimated.  The  equation 
involved  in  Ref.  6 was: 


where  the  primes  denote  differentiation  with  respect  to  a non-dimensional 
time,  T,  v(T)  is  non-dimensional  excitation,  and  the  non-dimensional 
"roll  angle"  6 is  the  ratio  of  roll  angle  (cp)  to  the  angle  ( cp  ) at  which 
Eq.  (6)  would  become  statically  unstable  if  £ is  negative;  that  is, 


Making  the  equivalent  non- 1 inear izat ion  of  Eq.  (b) , Eq.  (77)  is  approxi 
mated  by  Eq.  (6)  if  the  cp's  in  Eq.  (6)  are  replaced  by  0’s,  and: 


One  of  the  cases  computed  in  Ref.  6 was  selected  for  the  present. 
This  case  (No.  3210)  involved  values  of  linear  and  quadratic  damping 
(0,0)  thought  typical  for  a ship  with  bilge  keels.  The  particular  values 
of  damping  coefficients  were  a = 0.03  and  P = 1.0.  For  this  case  simu- 
lations were  carried  out  in  Ref.  6 for  three  levels  of  random  excitation. 
The  spectrum  involved  was  an  approximation  to  the  effective  wave  slope 
for  the  12th  ITTC  two  parameter  wave  spectrum.  In  Ref.  6 the  method  of 
increasing  level  of  excitation  was  to  multiply  a given  random  time  history 
by  a constant  factor.  Accordingly,  the  excitation  spectrum  for  each  level 
of  excitation  was  the  same  except  for  a constant  factor  in  the  spectral 
dens i ty. 
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Table  1 summarizes  some  of  the  parameters  involved  in  the  simula- 
tion as  well  as  the  present  computation.  The  nominal  rms  excitation  is 
indicated  for  the  three  levels,  as  is  the  rms  non-dimensional  roll  found 
in  the  simulation  for  each  excitation  level.  The  rms  roll  velocity  was 
estimated  from  the  observed  roll  spectra.  Because  the  roll  spectra  were 
narrow  band,  the  rms  roll  velocity  is  nearly  o times  the  rms  roll,  or,  in 
the  present  case  is  numerically  nearly  the  same  as  the  rms  ro'l.  Assuming 
on  the  basis  of  the  results  of  Ref.  6, that  the  roll  velocity  maxima  may 
be  expected  to  be  distributed  according  to  the  Rayleigh  distribution,  a 
reasonable  estimate  of  the  range  of  roll  velocity  (R)  was  throught  to  be 
2.5  times  rms  roll  velocity,  and  these  values  are  shown  in  Table  1.  (2.5 

rms  roll  velocities  corresponds  to  the  expected  value  of  the  1/10  highest 
ma  x i ma . ) 

Since  the  objective  of  this  first  evaluation  was  to  see  if  the 
formulae  developed  in  previous  sections  make  any  sort  of  sense,  no  itera- 
tions for  the  appropriate  value  of  the  parameter,  R , were  carried  out. 
Instead  the  value  estimated  from  the  known  simulation  results  (Table  1) 

was  substituted  in  Eqs  (78)  and  the  resulting  estimates  for  a and  a 

1 3 

are  given  in  the  table. 


TABLE  I 

Parameters  for  Case  3210  of  Ref.  6 


(or  = 0.03: 

q 

11 

(LI 

Excitation  Level 

1 

2 

3 

Excitation  rms 

0.004 

0.012 

0.036 

RMS  roll  found  from 
simulation 

0.018 

0.042 

0.09 

Estimated  value  of  R 

0.045 

0. 107 

0.236 

from  simulation 
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Eq.  (76)  was  programmed  for  evaluation,  given  the  coefficients  a, 

£,  or  , and  a and  a single-sided  excitation  spectrum  U ( u>)  defined  at 
13  x 

points  uniformly  spaced  along  the  frequency  axis.  The  linear  roll  spectrum 
required  by  Eq.  (76)  is: 

2 

U(u>)  = IG  (u>)  I U (w) 

1 * 

The  excitation  spectrum  used  in  each  evaluation  was  that  computed  from  the 

simulated  excitation  data  of  Ref.  6,  the  damping  coefficients,  a and  a 

1 3 

were  those  shown  in  Table  1. 

Figure  3 indicates  the  results  for  the  lowest  level  of  excitation. 

As  may  be  noted  in  the  Figure  the  excitation  spectrum  was  quite  broad 
relative  to  the  roll  response  (the  high  frequency  "tail"  is  omitted  from 
the  plot).  The  roll  spectrum  estimated  from  the  simulated  data  is  shown 
plotted  as  circles,  one  for  each  individual  estimate.  The  90$  confidence 
bounds  of  these  estimates  may  be  formed  by  adding  and  subtracting  about 
15$  of  each  estimate.  The  results  of  the  evaluation  of  Eq.  (76)  at  the 
same  frequency  interval  as  for  observation  are  shown  to  be  in  reasonably 
good  agreement  with  observation  — probably  well  within  the  statistical 
confidence  of  the  observations  since  in  addition  to  the  uncertainties  in 
the  observed  roll  spectrum,  there  are  uncertainties  of  the  same  magnitude 
in  the  excitation  spectrum. 

Also  indicated  in  Figure  3 is  U(  uu)  the  "linear"  spectrum,  and  the 
component  of  the  spectrum  resulting  from  the  first  of  the  three  terms  in 
Eq.  (76).  The  double  convolutions  in  Eq.  (76)  represent  a significantly 
greater  computation  effort  than  the  first  term.  It  may  be  noted  from  the 
figure  that  the  convolutions  evidently  cannot  be  neglected.  The  general 
relationship  between  the  linear  spectrum,  the  evaluation  of  Eq.  (7 6),  and 
the  contributions  of  the  convolutions  is  practically  the  same  as  that  shown 
by  Yamanouchi^.  An  evaluation  of  Eq.  (7 6)  was  also  made  with  £ = 0;  that 
is,  without  restoring  nonlinearity.  The  results  differed  insignificantly 
from  those  shown  in  Figure  3. 
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Figure  4 indicates  results  for  the  Level  2 excitation.  In  this 
case  the  rms  excitation  is  a factor  of  three  greater  than  that  of  Level  1, 
Figure  3*  The  deviations  from  observation  of  an  evaluation  of  Eq.  (76) 
for  this  case  are  stronger  than  those  shown  in  Figure  3,  but  possibly 
still  not  too  significant.  The  effect  of  removing  the  restoring  non- 
linearity are  also  shown  in  Figure  h.  Removal  of  this  nonlinearity 
emphasizes  a dip  in  the  spectrum  at  undamped  resonance.  The  dip  is 
caused  by  the  first  term  in  Eq.  (76),  the  contributions  of  the  double 
convolution  terms  do  not  make  up  the  deficiency. 

Figure  5 shows  the  results  for  Level  3 excitation,  (in  this  case 
rms  excitation  is  a factor  of  9 greater  than  that  of  Level  1,  Figure  3-) 
Here  the  results  of  evaluation  of  Eq.  (76)  differ  significantly  from  ob- 
servation. As  may  be  noted  from  the  figure,  the  results  of  omitting  the 
restoring  nonlinearity  come  close  to  observation.  It  appears  that  at  this 
level  of  excitation  something  is  missing  with  respect  to  the  computation 
of  the  effect  of  the  restoring  nonlinearity.  The  implication  is  that 
some  or  all  of  the  terms  resulting  from  the  fifth  degree  frequency  response 
(Eq.  73  through  75)  may  be  required.  This  last  (Level  3)  excitation  is 
relatively  severe.  The  observed  non-dimensional  rms  roll  (Table  l)  was 
roughly  9$  of  the  half-range  of  static  stability  of  Eq.  (77).  If  this 
magnitude  is  carried  over  into  ship  terms  as  in  Ref.  6,  rms  dimensional 
roll  might  correspond  to  6 or  7°,  and  if  the  Rayleigh  distribution  holds, 
the  average  of  the  10$  highest  out-to-out  rolls  would  correspond  to  30  or 
35° • This  is  usually  considered  severe  rolling,  but  it  is  not  by  any  means 
beyond  conditions  sometimes  experienced. 


It  appears  from  these  examples  that  the  development  for  the  spectrum 
does  make  sense  relative  to  the  results  of  the  simulations  of  Ref.  6.  The 
question  of  how  far  the  series  must  be  carried  out  so  as  to  be  valid  for 
the  entire  practical  range  of  ship  roll  has  not  been  answered,  but  it 
appears  possible  that  terms  involving  the  fifth  degree  response  function 
(Eq.  73  through  75)  may  have  to  be  retained  for  severe  rolling. 
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FIGURE  5-  COMPARISON  OF  OBSERVED  ROLL  SPECTRUM  FROM 
SIMULATION  WITH  RESULTS  OF  AN  EVALUATION  OF 
Eq.  76  FOR  LEVEL  3 EXCITATION 


On  the  other  hand,  the  percentage  errors  which  would  be  made  in 
the  prediction  of  rms  roll  according  to  Eq.  (7 6)  were  found  to  be  (rela- 
tive to  the  observed  value  from  the  simulation): 

Level  Figure  $ 

1 3 -11 

2 L - 7 

3 5 +18 

These  magnitudes  are  the  same  as  or  slightly  larger  in  magnitude  than  the 
90$  confidence  interval  on  rms  narrow  band  roll  which  would  be  implied  by 
the  sample  lengths  employed  in  the  Monte-Carlo  type  simulations  of  Ref.  6. 
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CONCLUDING  REMARKS 

The  point  of  undertaking  the  present  effort  was  to  see  if  the 
functional  series  model,  which  has  many  conceptual  attractions,  might  be 
applied  to  ship  rolling.  The  immediate  objective  of  the  work  was  to 
develop  means  of  estimating  the  scalar  spectrum  of  nonlinear  roll,  assum- 
ing that  the  usual  single  degree  of  freedom  equation  holds. 

In  order  to  apply  the  functional  series  to  the  usual  single-degree- 
of-freedom  equation  it  is  necessary  to  replace  the  quadratic  damping  term 
with  an  analytic  form.  Separately  reported  work  under  the  present  project 
has  indicated  that  it  is  realistic  to  replace  the  quadratic  term  by  a 
cubic  term.  Alternately,  it  appears  that  the  linear-plus-quadratic  damping 
of  the  usual  equation  may  be  approximated  by  a linear-plus-cubic  damping  — 
the  net  effect  of  this  "equivalent  non-linearization"  being  an  improvement 
over  simple  linearization  in  the  representation  of  linear-plus-quadratic 
damping  in  the  time  domain. 

Given  analytic  forms  for  all  the  non-linearities  in  the  equation, 
the  functional  series  may  be  applied,  and  the  derivation  of  the  associated 
series  of  frequency  response  functions  is  straightforward  if  laborious. 

Given  these  functions,  expressions  for  the  scalar  spectrum  of  roll  may  be 
derived  as  a series,  and  these  expressions  have  been  carried  through  terms 
corresponding  to  the  fifth  degree  in  the  basic  functional  series.  Owing 
to  the  nature  of  the  analytic  s ingle-degree-of-f reedom  equation,  the  terms 
of  even  degree  in  the  functional  series  were  found  to  be  zero.  Accordingly, 
the  representation  of  ship  rolling  and  the  expressions  for  the  spectrum 
involves  a functional  series  having  only  terms  of  odd  degree. 

A trial  evaluation  of  the  spectrum  of  non-linear  roll  was  made  for 
a case  which  had  previously  been  simulated  and  which  included  nonlinear 
restoring  moment.  Only  the  first  and  third  degree  terms  in  the  functional 
series  were  assumed  to  be  significant  in  this  evaluation.  The  results 
agreed  fairly  well  with  observation  for  moderate  to  low  level  rolling, 
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but  not  so  well  for  severe  rolling.  It  is  not  clear  from  those  evaluations 
which  could  be  completed  whether  or  not  it  will  be  necessary  to  include 
the  contributions  of  the  fifth  degree  response  functions.  The  errors  in 
estimation  of  rms  rolling  using  only  first  and  third  degree  terms  are  not 
impossibly  large. 

The  results  of  the  present  development  are  inherently  a series.  In 
practical  application  it  is  always  necessary  to  find  out  how  far  a series 
must  be  carried  in  order  to  achieve  a given  precision  of  estimate.  Though 
this  question  has  not  been  answered  by  the  present  study,  it  appears  pro- 
bable that  the  functional  series  for  rolling  need  be  carried  no  further 
than  terms  involving  fifth  degree  functionals,  and  possible  that  truncation 
after  the  third  degree  functional  may  suffice  in  most  applications.  If 
the  former  is  true,  the  complexity  of  the  model  may  well  be  more  than  can 
be  managed  in  practical  application.  Ho/ever,  if  the  latter  turns  out  to 
be  the  case,  and  non-linear  rolling  is  representable  by  a model  of  first 
and  third  degree  functionals,  practical  application  does  not  appear  so 
threatening  and  the  functional  series  approach  might  at  least  be  developed 
into  a viable  alternative  method. 


R-189U 


i 


I 


1 

\ 


r 


. 

t 


RECOMMENDATIONS 

It  should  perhaps  be  pointed  out  that  the  questions  about  the 
adequacy  of  a series  including  functionals  of  third  degree  relative  to 
a series  including  functionals  of  fifth  degree  are  themselves  all  relative 
to  the  assumption  that  beam  sea  ship  rolling  throughout  the  range  of  prac- 
tical interest  is  describable  by  a single  equation  which  involves  non- 
linearities  of  quite  specific  form.  This  same  assumption  has  been  made 
in  applying  alternate  approaches  to  the  problem  (statistical  linearization, 
perturbation  methods,  Fokker-Planck,  etc.).  It  seems  fair  to  speculate 
that  the  refinements  possible  in  any  of  the  approaches  may  over-reach  the 
validity  of  the  assumed  physical  model  when  rolling  throughout  the  range 
of  practical  interest  is  considered,  and  that  some  more  attention  paid  to 
the  fundamental  physical  nature  of  the  nonlinearities  and  their  coupling 
with  other  modes  of  motion  might  be  profitable. 

Considering  the  present  functional  series  approach  to  the  spectrum 
of  roll  as  defined  by  the  assumed  s i ngl e-degree-of-f reedom  equation,  many 
more  numerical  eval uations  wi 1 1 be  required  to  assess  the  range  of  validity 
and  the  highest  degree  of  functional  required  for  given  magnitudes  of 
nonlinear  coefficients  in  the  equation. 

Somewhere  in  the  development  of  alternate  practical  methods  of 
handling  nonlinear  rolling,  mul t i -degree-of-f reedom  systems  must  be  ac- 
commodated and  therefore  the  compatibility  of  the  functional  series 
approach  with  coupled  systems  will  need  to  be  investigated,  and  the 
implications  with  respect  to  experiment  be  clarified  so  that  realistic 
verification  experiments  can  be  designed. 
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roll  restoring  function 

restoring  coefficients,  Eq.  25 
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frequency  response  function  of  linear  output  filter 
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